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We outline the steps in a derivation of the statement that the SU(2) gauge theory is in a confining phase for all 
values of the coupling, < /3 < oo, defined at lattice spacing a. The approach employed is to obtain both upper and 
lower bounds for the partition function and the 'twisted' partition function in terms of approximate decimation 
transformations. The behavior of the exact quantities is thus constrained by that of the easily computable 
bounding decimations. 



1. Introduction 

A very large body of work has been performed 
by the lattice community in recent years in an ef- 
fort to isolate the types of configurations in the 
functional measure responsible for maintaining 
one confining phase for arbitrarily weak coupling 
in SU{N) gauge theories. Thick vortices disor- 
dering the vacuum over long scales at negligible 
local free energy cost have emerged as a primary 
mechanism (see ^ for a review and references). 
Nevertheless, a complete, direct derivation from 
first principles of this extraordinary and unique 
feature of SU{N) theories (shared only by non- 
abelian ferromagnetic spin systems in 2 dimen- 
sions) has remained elusive for three decades. 

The difficulty stems from the multi-scale na- 
ture of the problem: passage from short distance 
ordered perturbative regime to long distance dis- 
ordered non-perturbative confining regime. It 
can be addressed in principle only by a non- 
pertubative block-spinning procedure bridging 
short and long scales. Exact block-spinning 
schemes in gauge theories so far appear virtually 
intractable, both analytically and numerically. 

There are, however, approximate decimation 
procedures that can provide bounds on judicially 
chosen quantities. The idea is not new, but in 
the past only upper bounds were considered in 
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this context. The basic strategy in the follow- 
ing is to obtain both upper and lower bounds for 
the partition function and the partition function 
in the presence of a 'twist' (external center fiux). 
The bounds are in terms of approximate decima- 
tions of the 'potential moving' type, which can be 
explicitly computed to any accuracy. This leads 
to a rather simple construction constraining the 
behavior of the exact partition functions, and, 
through them, the exact vortex free energy and 
other order parameters, by that of the bounds. 
They thus are shown to exhibit confining behavior 
for all values of the inverse coupling, < (3 < oo, 
defined at lattice spacing a (UV cutoff) held fixed. 
Only the SU{2) case is considered explicitly here, 
but the same development can be extended to 
general SU{N). 

2. Decimations 

We begin by some plaquette action, e.g the Wil- 
son action Ap{U,a) — ^ KetiUp, defined at lat- 
tice spacing a. The character expansion of the 
exponential of the action: 

F{U, a) = e^-(^) = P^^P^ «) '^^ ^^(f^) (1) 
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is given in terms of the Fourier coefficients: 

Fj = J dU F{U,a)^xdU) ■ (2) 

Here Xj denotes the character of the j-th rep- 
resentation of dimension dj. Thus, for SU(2), 
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0, i 1, S, . . ., and dj — {2j + 1). In terms 



of normalized coefficients: 

one then has 
FiU,a) = Fo 



(3) 



i^o/(C/,a) 



(4) 



For a reflection positive action one necessarily 
has: 

Fj >0 , hence 1 > > , aU j .(5) 

The partition function on lattice A is then 



Za{P) = F, 



A| 



dUA 



(6) 



We now consider RG decimation transforma- 
tions a Xa. This involves partitioning the 
lattice in c?-dimensional decimation cells of side 
length Xa. Simple approximate transformations 
of the 'potential moving' type are implemented 
by 'weakening', i.e. decreasing the Cj's of pla- 
quettes interior to the cells, and 'strengthening', 
i.e. increasing Cj's of cell boundary plaquettes. 
The simplest scheme [Sj, which is adopted in the 
following, implements complete removal, Cj — 0, 
of interior plaquettes. This may be pictured as 
moving the interior plaquette interactions to the 
cell boundaries. The operation can be decom- 
posed into elementary moving steps along each 
positive direction as illustrated for, say, the x^- 
direction in Figure ^ The (A — 1) interior plaque- 
ttes (shaded) are moved and merged with the cor- 
responding boundary plaquette (bold) into one 
boundary plaquette with renormalized interac- 
tion Ap{u) CMu)- 

This basic operation is successively performed 
in all directions. Note that, in d dimensions, there 
are {d — 2) normal directions into which a plaque- 
tte can move. A plaquette moved to the {d — 1)- 
dimensional cell boundary can still be moved in 
{d — 3) directions inside the cell boundary. The 
moving operation terminates when all plaquettes 
have been moved in this manner to form a tiling 
of the 2-dimensional faces of a lattice of spacing 




X = 2 



Figure 1. Basic plaquette moving operation. 



Xa. The integrations over the bonds of the tiling 
plaquettes inside each such face can now be per- 
formed exactly (being 2-dimensional) , thus merg- 
ing the tiling plaquettes into one plaquette of the 
coarse lattice of spacing Xa. 

In terms of the definitions and notations in- 
troduced above, the end result can be concisely 
stated as follows. Under successive decimations 



a - 
A 



Xa X^a 

. Ad) ^ A(2) 



^ a(" 



the resulting RG transformation rule is: 

/(f/,n-l)->Fo(n)/([/,n) 

with 



fiU,n) 



and 

Cj{n) = Cj(n)^' , 
where 

Cj{n)=F,{n)/Fo{n) < 1 



Fo{n)^Fo{ny 



dU 



(7) 
(8) 

(9) 

(10) 
(11) 



The renormalization parameter controls by how 
much the plaquettes remaining after each decima- 
tion step have been strengthened to compensate 
for the removed plaquettes. What has been con- 
sidered in the literature almost exclusively is the 
choice C = A'^'^^^'. This is essentially the original 
choice in P], and will be referred to as the MK 



3 



choice. Here we consider C an adjustable param- 
eter. 

The resulting partition function after n deci- 
mation steps is: 



|A|/A" 



rfC/AW n MU,n). (12) 

As a point of notation, the dependence of the 
quantities _fo(rt), Cj{n), etc. on variables such as 
A, ^, or the set of couplings /3, which characterize 
the choice of decimation and action at the original 
spacing a, will not be indicated explicitly unless 
specific reference to it is required. 

It is important to note that after each decima- 
tion step the resulting action retains the original 
one-plaquette form but will, in general, contain 
all representations: 

Ap{U,n)^boin) + Y, P) X^U) ■ (13) 

Furthermore, among the effective couplings 
/3j(n, (3) some negative ones may in general occur. 
These features are present even after a single deci- 
mation step a — *■ Aa starting with the usual single 
representation (fundamental) Wilson action. 

Preservation of the one-plaquette form of the 
action is of course what makes these decimations 
simple to explore. The rule specified by (TJl- (fTT|l 
is meaningful for any real (positive) Here, how- 
ever, a basic distinction can be made. For inte- 
ger ^, the important property of positivity of the 
Fourier coefficients in (Q, 10}: 

Fo(7i)>0 , c,(n)>0, (14) 

and hence refiection positivity are maintained at 
each decimation step. This, in general, is not the 
case for non-integer Thus non-integer C results 
in approximate RG transformations that violate 
the reflection positivity of the theory (assuming 
a reflection positive starting action).'^ From now 
on we assume that (|14ll holds. 

^It is worth noting in this context that various numerical 
investigations of the standard MK recursions, at least for 
gauge theories, appear to have been carried out, for the 
most part, for fractional A, (1 < A < 2), which corresponds 
to non-integer e.g. see UJ. 



There are various other interesting features of 
such decimations. The following property, in par- 
ticular, is important. One can show that, given 
the coefficients Cj{n) after n decimations, one has 
the general relation 



(15) 



It follows from 115|) that the {I2) norm of the vec- 
tor c(n-|- 1) formed from the Cj{n-\-l) coefficients 
is bigger than that of the vector of the Cj (n) : 



||e(n + l)||2>||c(n)||2 



(16) 



All coefficients being positive, this implies that 
(|15ll must hold also component-wise, ie. 



Cj{n + 1) >Cj{n) , 



(17) 



for at least a subset of components giving the 
dominant contribution to the norms. For ^ = 1, 
it follows immediately from ifTUI) - lfTT|) that lfT7|l 
holds as an equality for all j. For general one 
finds by explicit numerical evaluations that, in 
fact, H17|) also holds for all j. This can be shown 
analytically in special cases and can probably be 
proved in general. 

Alternatively, instead of attempting such a 
proof, one may impose (|17|l as part of the spec- 
ification of the decimation transformation. This 
amounts to replacing 110() by 



Fo{n) LFo(n) 



+ Cj(n- 1)9 



Cj{n-1) 

Fo{n) 



(18) 



3. The exact partition function 



Since our decimations are not exact decimation 
transformations, the partition function does not 
in general remain invariant under them. The sub- 
sequent development hinges on the following two 
basic propositions that can now be proved: 

(I) With C = X'^^^: 



ZAiP.n) < ZA(/3,n + l) 
(II) With C = 1: 

ZA{(3,n + l) < ZA{(3,n) 



(19) 



(20) 
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Note that for d = 2 - (EOl express the well- 
known fact that the decimations become exact. 
For d > 2, in both (I), (II) one in fact has strict 
inequality. 

(I) says that modifying the couplings of the 
remaining plaquettes after decimation by taking 
( = A'^"^ (standard MK choice [2]) results into 
overcompensation (upper bound on the partition 
function). (II) says that decimating plaquettes 
while leaving the couplings of the remaining pla- 
quettes unaffected results in a lower bound on the 
partition function. Translation invariance, con- 
vexity of the free energy and reflection positivity 
underlie fH^ . 

Consider now the, say, (n — l)-th decimation 
step with Fourier coefhcients Cj{n — 1), which we 
relabel Cj(n — 1) = Cj{n — 1). Given these Cj{n — 
1), we proceed to compute the coefficients Fo(n), 
Cj{n) of the next decimation step according to 
©-(HU above with C = A''"^ 

Then introducing a parameter a, (0 < a), de- 
fine the interpolating coefficients: 

c,(n,a) = c;,(n-l)^'(i-")c,(n)", (21) 

so that 

Cj{n) : a — I 

Cj{n~l)^^ : a = 

The a = value is that of the n-th step coeffi 
cients resulting from (P|l- (|ll|l with C = 1- 

Thus defining the corresponding partition func 
tion 



Cj(n, a) 



n-l 



A|/A" 



m=0 



J dUj^ir., [] /p(C/,n, a) (23) 



where 
fp{U,n,a) 



(24) 



we have from (Unj, ^^), and above: 

Za(/3, 0, n) < Za(/3, n - 1) < Z^W, 1, n) . (25) 

Now the partition function (|23|l is a continuous in 
a. So H25(l implies that, by continuity, there exist 
a value of a: 

a = «(")(/?, A, A), < a(")(/3,A,A) < 1 



such that 

ZA(/3,a("),n) = ZA(/3,n-l) 



(26) 



In other words there is an a at which the n-th 
decimation step partition function equals that ob- 
tained at the previous decimation step; the parti- 
tion function does not change its value under the 



decimation step A" 



A"a. 



So starting at original spacing a, at every dec- 
imation step TO, (m = 0, 1, ■ • ■ ,n), there exist a 
value < a^'") < 1 such that 



Za(/?, ^("+1), m + 1) = Za(/3, m) . 



(27) 



This then gives, after n successive decimations, 
a representation of the exact partition function in 
the form: 

Za(/3) = fI^I J dUA n fp{U,a) 

n 



pGA 
v'^'lAI/A" 



ZAw(/3,n,a(")) 



(28) 



where 



(22) ZA(„,(/3,n,a(")) 



dU, 



A(") 



n /p(C/,n,a(")),(29) 

peA(") 



i.e. a representation in terms of the successive 
bulk free energy contributions from the a — > A — > 
• • • — > A"a decimations and a one-plaquette effec- 
tive action on the resulting lattice A^"). 

At weak and strong coupling a'™-* may be es- 
timated analytically. At large /3, where the dec- 
imations approximate the free energy rather ac- 
curately, the appropriate a values are very close 
to unity. At strong coupling they may be esti- 
mated by comparison with the strong coupling 
expansion. On any finite lattice there is also a 
weak volume dependence as a correction which 
goes away as an inverse power of the lattice size. 

For most purposes the exact values of the 
Q,(m)ig^ beyond the fact that are fixed between 
and 1, are not immediately relevant. The main 
point of the representation (|28|l is that it can in 
principle relate the behavior of the exact theory 
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to that of the easily computable approximate dec- 
imations. 

Indeed, starting from the Cj{n — 1, a'"^^-')'s at 
the (n — l)-th step, consider the coefficients at 
the next step, and compare those evaluated at 
a = Q!^"-*, i.e. Cj{n, a = a'"-*), to those evaluated 
at a = 1, i.e Cj{n,a = 1) = Cj{n). The latter 
will be referred to as the MK coefhcients. (Recall 
that a = 1 C = the standard MK 

choice. The absence of a tilde on a coefficient in 
the following always means that it is computed 
at a — 1.) According to (I), the MK coefficients 
give an upper bound. 

Now, ((SH) implies that 

Cj{n,a) < (l-a)cj(n-l)^' +aCj{n), (30) 

from which, by property H17|) . one has 

Cj(n, a) < Cj{n) for any < a < 1 . (31) 

This has the following important consequence. 

H31() says that the Fourier coefhcients of the 
representation H28I) are bounded from above 
by the MK coefhcients (a = 1). Thus, 
if the Cj(n)'s are non-increasing, so are the 
Cj{n, a). The Cj (rt)'s must then approach a fixed 
point, and hence so must the Cj(n, a)'s, since 
Cj{n),Cj{n,a) > 0. Note the fact that this con- 
clusion is independent of the specific value of the 
a's at every decimation step. 

In particular, if the Cj(ri)'s approach the strong 
coupling fixed point, i.e. Fq ^ I, Cjin) — > 
as n — > oo, so must the Cj(n, q;)'s of the exact 
representation. // the MK decimation coefficients 
flow to the confining regime, so do those in the 
exact representation 1^2 1^) . As it is well-known 
by explicit numerical evaluation, the MK deci- 
mations for SU{2) and SU{3) indeed flow to the 
strong coupling confining regime for all /? < oo 
and d < 4. Above the critical dimension d — 4, 
the decimations result in free spin wave behavior. 

What do these results imply about the ques- 
tion of confinement in the exact theory? Though 
strongly suggestive, the fact that the long dis- 
tance part, Z/^(„), in 128|) flows in the conflning 
regime does not suffice to answer the question. It 
is the combined contributions from all scales in 
(|28|l that combine to give one bulk quantity, the 



exact free energy — In Z\{P). As it is well-known, 
it is not possible to unambiguously determine the 
long distance behavior of the theory from that of 
a bulk quantity like the free energy. For that one 
needs to consider appropriate long distance order 
parameters. 

4. Order parameters - Vortex free energy 

The above derivation leading to the represen- 
tation H28|l for the partition function cannot be 
applied in the presence of observables without 
modiflcation. Thus, in the presence of operators 
involving external sources, such as the Wilson or 
't Hooft loop, translation invariance is lost. Re- 
flection positivity is also reduced to hold only in 
the plane bisecting the loop. Fortunately, there 
are other order parameters that can characterize 
the possible phases of the theory while maintain- 
ing translational invariance. They are the well- 
known vortex free energy, and its Z{N) Fourier 
transform (electric flux free energy). They are in 
fact the natural order parameters in the present 
context since they are constructed out of parti- 
tion functions. Recall that the vortex free energy 
is defined by 

g-F„(r) ^ z^{t)/Za . (32) 

Here Z\{t) denotes the partition function with 
action modified by the 'twist' r G Z{N) for ev- 
ery plaquette on a coclosed set of plaquettes V 
winding through the periodic lattice in {d ~ 2) 
directions, i.e. through every [//j/]-plane for fixed 

Ap{Up) ^ ApirUp) , if p e V . 

One has Z\{1) = Z\. A nontrivial twist {t ^ 1) 
represents a discontinuous gauge transformation 
on the set V which introduces Tri{SU{N)/Z{N)) 
vortex flux rendered topologically stable by be- 
ing wrapped around the lattice torus. As indi- 
cated by the notation, Z\(t) depends only on 
the presence of the flux, and is invariant under 
changes in the exact location of V. The vortex 
free energy is then the ratio of the partition func- 
tion in the presence of this external flux to the 
partition function in the absence of the flux (the 
latter is what was considered above). As it is 
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well-known the possible phases of a gauge theory 
(Higgs, Coulomb, confining) can be characterized 
by the behavior of (|32|l . Furthermore, by rigor- 
ous correlation inequalities H^, the Wilson loop, 
Wilson line correlators, and the 't Hooft loop can 
in turn be bounded by the vortex free energy and 
its Z{N) Fourier transform. For SU{2) the only 
nontrivial element is r = — 1. 

The above development, in particular the 
derivation of (|28l) . should then be repeated for 
Za{—1) = Zj^ \ There is, however, a technical 
complication in obtaining the analog to (|28(l for 
Z'^\ The presence of the flux reduces reflection 
positivity to hold only in planes perpendicular to 
the directions in which the flux winds through 
the lattice. This may be circumvented by con- 
sidering Z\ + Zj^^ instead of Zj^\ Indeed, for 
Z\+Zj^'^'^ reflection positivity is easily checked to 
again hold in all planes. The proof of propositions 
(I) and (II), and the subsequent derivation lead- 
ing to (|2ZIl, and ^ can then be carried through 
for the quantity Za + Zj^ ' to obtain 

n 

Zdfi) + Z^-\fi) = n i^oM"""'!^!/^'"'' 

m=0 

• ( (/3, n, ) -f 47i, (/?, n, «("))) (33) 
where 

peA(") 

■ n /p((-l)C/,^,«'"^)). (34) 

pev 
peA(") 

The values {a^™^} in H33|l . as flxed by the inter- 
polating argument between the upper and lower 
bounds for the quantity Za -|- \ are a priori 
distinct from those in H28() flxed by the analogous 
argument for the quantity Z;^. It is not hard to 
see, however, that in fact they have to coincide 
for large lattices. 

Using I^Hl) and (j^Sl in (Ell then gives for the 
vortex free energy: 



It manifestly couples only to the long distance 
part as it should. Bulk (local) free energy con- 
tributions resulting from each successive decima- 
tion are unaffected by the presence of the flux 
and cancel. By our previous considerations, the 
strategycoefficients Cj(n, a*-"-') occurring in Zj^(n) , 
Z^^i) i'^ are bounded (regardless of the spe- 
ciflc values of the a^"^ parameters) by the corre- 
sponding MK coefflcients Cj(n) = Cj{n,a — 1). 
By taking n large enough, the latter tend to zero 
no matter how large the initial /3 is chosen. Thus 
taking n large enough to enter the strong coupling 
regime, (|35f) can be evaluated exactly within the 
convergent strong coupling cluster expansion in 
terms of the Cj(n, a("-))'s by standard computa- 
tions. The vortex free energy is thus explicitly 
shown to exhibit confining behavior. This, by 
the known inequalities Q relating the vortex free 
energy to the Wilson loop, then also necessarily 
implies area law for the latter. 

An expanded account with the proofs of (I), 
(II) and various other statements above will ap- 
pear elsewhere. 

I am grateful to the organizers of the QCD 
Down Under Workshop for organizing such a 
stimulating meeting in a thoroughly enjoyable en- 
vironment in Adelaide and the Barossa Valley, 
and to the participants for a great many discus- 
sions. 
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